PROPERTIES OF LINEAR INTEGRAL EQUATIONS RELATED TO 
THE SIX- VERTEX MODEL WITH DISORDER PARAMETER II 
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Abstract. We study certain functions arising in the context of the calculation 
of correlation functions of the XXZ spin chain and of integrable field theories 
related with various scaling limits of the underlying six-vertex model. We show 
that several of these functions that are related to linear integral equations can be 
obtained by acting with (deformed) difference operators on a master function <J>. 
The latter is defined in terms of a functional equation and of its asymptotic behavior. 
Concentrating on the so-called temperature case we show that these conditions 
uniquely determine the high-temperature series expansions of the master function. 
This provides an efficient calculation scheme for the high-temperature expansions 
of the derived functions as well. 



1. Introduction 

In [4,6,7, 10] an algebraic structure of the correlation functions of the XXZ model 
(or of the underlying six- vertex model) called factorization was identified in a rather 
general setting. It was shown that the correlation functions consist of an algebraic 
part and of a physical part. The algebraic part is determined by a set of operators 
t*, b, c, b*, c* which act on a space of quasi-local operators. They do not 

depend on the physical parameters like temperature, magnetic field etc. The operators 
b, c (annihilation operators) and b*, c* (creation operators) are two pairs of Fermi 
operators satisfying characteristic anti-commutation relations. The creation operator 
t* is 'bosonic' and commutes with all Fermi operators. With the help of the creation 
operators one can construct a 'fermionic basis' of the space of quasi-local operators. 
In contrast to the algebraic part the physical part depends on the physical parameters 
and is represented by two transcendental functions p and w. 

The proof of the factorization in [10] covers the rather general case of a six-vertex 
model which is finite in vertical, so-called Matsubara direction and carries arbitrary 
inhomogeneity parameters on the horizontal lines. The proof is based on the properties 
of the fermionic and bosonic operators. It also relies on a generalization of the density 
matrix (denoted Z K in [10]), such as to include a disorder parameter a, and on the 
use of so-called (/-deformed Abelian integrals of the second kind. The transcendental 
functions p and co which constitute the physical part also depend on one or two 
spectral parameters (,C- The function p(() is the ratio of two maximal transfer 
matrix eigenvalues with different twist parameters. It is related to the one-point 
correlator. The function u((, £') is related to the nearest-neighbor correlators, but its 
definition is more involved. The authors of [10] demonstrated that it is determined by 
a 'normalization condition' for their (/-deformed Abelian integrals of the second kind. 

In [2] we suggested an alternative description of the function uj by means of the 
solutions of certain linear and non-linear integral equations. We also proved the 
equivalence of this description and the one suggested in [10]. In [3] we discussed 
some properties of those integral equations in the case of finite temperature. The 
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normalization condition for oj turned out to be a consequence of the integral equations 
for the auxiliary functions a and G introduced in [2,8, 15, 16]. An important point 
in [3] was to introduce the 'dressed charge' a which, like the function G, fulfills a 
linear integral equation. Using Baxter's TQ-relation, it appeared to be possible to 
find a solution to this equation in 'explicit form' in terms of the ratio $ of eigenvalues 
of two Q-operators with different twist parameters. 

In our recent works [1,5] we studied the scaling limit towards CFT on a cylinder. 
We conjectured that the creation operators t*, b*, c* have well-defined scaling limits 
r*, f3* , 7*. We could identify the states generated by their action on the vacuum with 
the Virasoro module of the CFT descendant states up to level 8. This identification 
was achieved by considering the three-point CFT correlators on the cylinder with a 
descendant field P o ,({l_ n })^ Q (0) at the origin and two asymptotic 'edge'-fields (f> K (oo) 
and 0_ K /(— oo) at ±oo or their descendants, where the primary field 4> a has conformal 
dimension A a = v 2 a(a — 2)/(4(l — v)). On the other hand we took the scaling limit 
of the determinant formula derived in [10] for the generating function of correlation 
functions of the XXZ model, which is realized as a six-vertex model on a lattice 
with cylindrical topology and with creation operators inserted. This involves, in 
particular, taking the scaling limits p sc and u sc of p and uj. Comparing the coefficients 
of the asymptotic expansions at large spectral parameter, we could perform the above 
identification. Unfortunately, this was possible only in a weak sense, i.e. modulo 
integrals of motion. 

Let us further comment on the latter point, because it is of particular importance 
for the motivation of the present paper. First, we considered the six- vertex model at 
finite disorder parameter a and magnetic field k which also plays the role of a twist or 
flux parameter. In fact, in order to determine the correlation functions on the lattice 
in such a way that they have appropriate scaling limits, we also needed to introduce a 
kind of 'lattice screening operators' that change the spin by some number s. Then we 
came to a picture in which the twist parameter on the 'left' edge of the cylinder is k 
and on the 'right' edge k' = k + a + 2s(l — v)/u, where v is related to the deformation 
parameter q = e vv% . We argued that k' becomes an independent variable in the scaling 
limit. Thus, expectation values of quasi-local operators turn into three-point functions 
of two independent primary fields at ±oo and a descendent field at as was pointed 
out above. When k! = k it follows that p = p sc = 1, and the integrals of motion do 
not contribute. In this case we managed to develop (see [1,5]) a technique for the 
calculation of the coefficients in the asymptotic expansion of the function uj sc based 
on Wiener-Hopf factorization. 

The problem is how to extend this to the general case n' ^ k and p ^ 1. It seems 
we have to develop a new technique here. This work is a first attempt in this direction. 
We consider the temperature case, where the high-temperature expansion is a powerful 
tool to generate explicit results. 

We shall take up some ideas of the paper [3], where we studied the properties 
of linear integral equations, like the one for the dressed charge a, using the above 
mentioned function <3? in the temperature case. The idea elaborated below is to 
introduce a generalized function <I> of two arguments in such a way that it satisfies a 
functional relation which involves the original function <I> of one argument. Inserting 
the high-temperature expansions of the two functions we see that the expansion 
coefficients of the new function are recursively determined by those of the old function. 
This provides an efficient scheme to calculate them on a computer. 

We also consider the calculation schemes for the coefficients in the high-temperature 
expansions of the resolvent R, the function G, its dual G and the function of our main 
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interest oj. We find that they are all connected to our new function $ in a simple way: 

G(C, = A C F C ,$(C, CO , G(C, CO = - 2^fl- c A f ^(C, CO , 

where A and are certain difference operators defined in (2.10), (2.31). These 
formulae are valid in general, in particular, in the temperature case and in the scaling 
limit. This is the main result of this paper. The function 5>(C, CO has simpler analytic 
properties than the resolvent or the w-function and together with the function p 
contains the whole information that we need for the correlation functions. In this 
paper we show how to calculate the high-temperature expansion of $(C, CO- F° r now 
this seems to be easier as compared to the problem of calculating the asymptotic 
expansions in the scaling limit. We do not need to use the Wiener-Hopf technique 
here. Still, we plan to obtain the asymptotic expansion of the <E>-function in the scaling 
limit in our future work. Although our primary aim are CFTs this may turn out to 
be useful in conjunction with the recent progress in the understanding of the massive 
integrable Sine-Gordon theory [11-14] as well. 

The paper is organized as follows. In Section 2 we recall some definitions from the 
papers [2,3,5] and introduce those objects that we need for the further considerations. 
We also discuss the basic thermodynamic functions and the integral equations they 
fulfill. In Section 3 we define the function <3?(C) of one spectral parameter and recall 
how the 'dressed charge' a can be expressed in terms of this function. We derive an 
equation for the coefficients of the high-temperature expansion of 3>(C) which can 
be solved by iteration. In Section 4 we generalize this method to solve the linear 
equation for the resolvent R((, (')■ To this end we introduce a function CO which 
is a generalization of 3>(C)- We constitute Lemma 4.1 on the representation of the 
resolvent in terms of the function <&(C, CO an d obtain its high-temperature expansion. 
In Section 5 we formulate Lemma 5.1 comprising the above formulae representing the 
functions G, G and ui by the action of difference operators on <3?(C, CO- We also argue 
that these formulae are valid in general and not only in the temperature case. In the 
Appendices we prove Lemma 4.1 and 5.1 and show several lowest order coefficients of 
the high-temperature expansions of the above functions. 



2. Basic objects and equations 

2.1. Correlation functions and the functional Z K,S . This work is about the 
properties of special functions arising in the context of the calculation of correlation 
functions of the XXZ chain and of conformal field theories. Before defining these 
functions we would like to briefly sketch the context. For more details the reader is 
referred to [5,7,10]. 

A correlation function is an expectation value of a local operator 0, typically 
calculated as a thermal average by means of the statistical operator of, say, the 
canonical ensemble. Having in mind the scaling limits towards conformal and massive 
field theories it is natural and useful to generalize both, the notion of a local operator 
and the notion of the statistical operator. Instead of local operators so-called quasi- 
local operators with tail were introduced in [6,7]. In these articles the XXZ chain, a 
model of locally interacting spins, was considered on an infinite lattice. On such a 
lattice the action of S(k) = \ Y^j=-oo a j-> where a z is a Pauli matrix, makes sense. A 
quasi- local operator is an operator of the form q 2aS ^0, where is local. Expectation 
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values of local operators (of spin zero) can be defined by means of a functional 

Tr s Tr M (^ s) rs,M <? 2k5(oo) b* : • • • b* {q 2a 



(2.1) z K ' s i q 2aS Wo I = L : : ^ 

1 J TrgTrMl^^Ts.M g 2K ^ ■ ■ ■ b^ )0 (^(o)) } 

generalizing the canonical ensemble average of statistical mechanics. The 'lattice 
screening operators' b^j, which are the coefficients in the expansion of the singular 
part of b*(C) at C 2 = 0,* increase the spin of by s. This is compensated by the 
operator at the boundary which carries spin — s and ensures that the ice-rule of 

the six- vertex model is satisfied. As was discussed in [5], the functional (2.1) does not 

(—s\ 

depend on the concrete choice of ' . The spin is needed as an additional parameter, 
necessary to exhaust the full space of descendents in the conformal limit. 

The 'thermodynamic properties' of the average are determined by the monodromy 
matrix Ts,m and the operator q 2KS (°°) . Mathematically the monodromy matrix is 
defined by evaluating the universal iZ-matrix of t/^sfe) on the tensor product of two 
evaluation representations $]$ and .f)M- Both of them realize the space of states of 
a spin chain. With i^s we associate an infinite chain and with f)M a finite chain of 
length N, 

oo N 
Ts,m Yi T i> M ' T 3< M = t jm(1) ' t jm(0 Yl ^j,m(C/£m) ■ 

j=—oo m=l 

Here the £j € C are inhomogeneity parameters, and L is the standard L-operator of 
the six vertex model 

wo = q-^ aL - cM 9 * 4 - c(« - + ■ 

Different physical realizations of the spin chain (e.g. finite temperature or finite 
length) can be realized by appropriate choices of the parameters £j and k [2]. These 
parameters are also at our disposal for realizing scaling limits towards CFT [5] or 
massive integrable quantum field theories [12-14]. 

For the direction of infinite extension of the lattice one can change the boundary 
conditions and, instead of taking the traces at the right hand side of (2.1), insert 
two one-dimensional projectors \k)(k\ and \k + a — s, s){k + a — s, s\ at the boundary. 
Then 

r2 2 ^ 7K , s f«(0)J _> (« + "-', s\T SM <^ S • • • b^ Q ( g ^(0)O) \ K ) 

{ -> T J ( K + a - s , s |r S)M ^b^_ 1 ---b^ (^5(0))| K) » 

where \k) is the eigenvector of the transfer matrix Tm((,k) = Trj (Tj^mQ^ 3 ) with 
maximal eigenvalue T(£, At) in the spin-zero sector, and the eigenvector \k + a — s, s) 
corresponds to the maximal eigenvalue T(£, k + a — s, s) of the transfer matrix 
^m(C) k + a — s) in the sector with spin In [5] we argued that the combination 

(2.3) k = K + a + 2^s 

becomes an independent parameter in the CFT scaling limit. 



*For the case s < one can replace the operators b^, j by c^ j as in [14]. 

^For s 0, depending on the choice of the inhomogenieties, the eigenvalues may be generally 
degenerate. This is, however, irrelevant for our discussion below. 
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2.2. Functions p, uj and basic thermodynamic functions. In [7] the 'fermionic 
basis' described in the introduction was constructed. The creation operators t*, b* 
and c* generate the space of quasi-local operators by acting on the pure tail q 2aS (°) 
which plays the role of the Fock vacuum for these operators. Note that the creation 
operators are very special. Their most important property is their compatibility with 
the above defined functional Z K ' S revealing itself in the following formulae proved 
in [10] 



(2.4) Z^{t*(0(X)} = 2p((\ K , K + a,s)Z^{X}, 

Z«' s {c%0(X)} = -^Ju(tX\K,cx,s)Z"> s {b(0(X)}- 



Z^{b*(C)(X)} = u((,t\K,a, S )Z^{c(0(X)} de 



r 

d_e 

r 

Here the contour T encircles all the singularities of the integrand except £, 2 = ( 2 . The 
functions p and uj appearing in this theorem are in the center of interest of this work. 
We shall provide a precise definition below. They are fundamental for the description 
of all static correlation functions of the XXZ chain and its various scaling limits, since 
Wick's theorem combined with (2.4) implies the determinant formula 

(2.5) Z">'{t*«2) ■ ■ ■ t%>*(Ci + ) • ■ -b*(C + )c*(C-) • ■■c*(Ci)(q 2am )} 

v 

= T[M(i\ K , K + a , s ) x det f^(C + ,Ck,a,s)) 

. "T \ J /i,j=l,--,r 

1=1 

From this formula one can obtain any correlation function by Taylor expanding both 
sides and comparing coefficients, because the operators appearing on the left hand 
side generate a basis of the space of quasidocal operators as was proved in [4] . 
The function p in (2.4) is the ratio of two eigenvalues of the transfer matrix 

foa\ (t\ , A T((,K + a- s,s) 

(2.6) p((\K,K + a,s) = — . 

In the following we shall replace it with 

T(O') 



(2.7) p(Ch»') = 



T(C,«) 



and treat k' as real parameter independent of a (which reappears below in the 
definition of uj). As was explained in Section 4 of [5] this replacement is possible and 
necessary in the CFT scaling limit. 

For the function uj the following representation was originally obtained* in [2] and 
then used in [5], 

(2-8) -w(C, C'l«, «) = (/left * /right + /left * R * /right) (C + MC, C» • 

In the CFT scaling limit the function uj((, ('\k, k'; a) is identical to w(C, C'I K > a > s ) m 
(2.4) if (2.3) is fulfilled [5]. 



*In [2] we used a slightly different notation. In particular, up to some multiplier, the first two 
terms in brackets on the right hand side of (2.8) were denoted 
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The definition of the functions entering into (2.8) is slightly involved. But they 
either derive from the elementary function 

(2-9) ^ a) = C ^^l) 

or from the functions appearing in the TQ-relation. For their definition we shall also 
need the two difference operators 

(2.10) A C /(C) = f(Cq) ~ fiCq' 1 ) , a c "/(0 = /(C?) " P(C\K, «')/(C) • 
Then 

(2.11) / left (C, C', «) = ^^(C/C, «) . /right (C, C, «) = <^(C/C, a) 
and 

(2.12) Wo (C, Cl«) = A" V(C/C, «) . 
where A -1 ^ is defined as a principal value integral 



(2.13) Ajck.n) /,r / -.- r > 0. -- < Ren < 0. 



d?? 2 V>(7y,a) ( 2 > o -- 

o 2vrir/2 2l /(l + (C/»7)-) ' ' " 

The function R is the resolvent of a linear integral operator. For its definition we 
introduce the kernels 

(2-14) K a (0 = ^.A^(C,a), K(0 = K (C). 

Then R is defined by the integral equation 

(2.15) R- R*K a = K a . 

Here and in (2.8) we used the notation 

(/*0)(C,C) = / dmir)) f((,r,)g( V X') 

7 

for the convolution. The 'measure' dm is given as 

dry 2 



(2-16) <!'»('/) , ; ,w, , , 

and the contour 7 is described below. 
The auxiliary function a is defined by 

(2 - 17) a(C ' K) = a(CMrt,-) 

where Q((, k) is the eigenvalue of the Q-operator (for the definition of the Q-operator 
see Section 3 of [5]). It satisfies Baxter's TQ-equation 

(2.18) T(C, k)Q(C, «) = d(C)Q(9C, «) + aiOQiq-'C, «) 
where g = exp (70^) and 

JV A 

(2.19) a(C) = n( 1 -9CV^ 2 ), d(C) = n( 1 "rt 2 /eI 



For our convenience we assume TV to be even. The function C K Q(C> K ) is a polynomial 
and, hence, is determined by its zeros Q (the Bethe roots) satisfying the Bethe ansatz 
equation 

o(C-,«) = -l, j = !,-■■, N/2. 
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Instead of (2.17) one can use a non-linear integral equation to characterize the aux- 
iliary function a [15, 16]. Its precise form depends on the position of the inhomogeneity 
parameters [2]. The non-linear integral equation is particularly useful for numerical 
calculations, for the calculation of the large N asymptotics and for performing the 
limit N — > oo in the temperature case 

(2.20) 6,-1 = exp (ttui/2 - p/N) , £ 2j = exp (-irui/2 + p/N) , j = l,...,N/2. 

In [5] we considered the conformal limit. In the present paper we deal with the 
temperature case with the inverse temperature (3 = T~ l . Then, after taking the limit 
N — > oo, the non- linear equation for the auxiliary function becomes 

(2.21) log(c(C, «)) = -2ivn/K - i sin {w)0e(Q - j ^ K((/0 log(l + a(£, «)) 

7 

with the 'bare energy' e(() = e(exp A) = coth (A) — coth (A + wvi). The contour 7 in 
the complex £ 2 plane encompasses the essential singularity at £ 2 = 1, but the poles of 
the kernel at £ = (q^ 1 are outside. In terms of the variable fi = log£ it is depicted 
in [2] with the opposite integration direction. In the following we assume k and v 
to be real. This corresponds to the unphysical case of a purely imaginary magnetic 
field, but is convenient for our purpose of studying the function uj by means of the 
high-temperature expansion. 

2.3. Functions G and G. Originally the correlation functions in the temperature 
case for a = were represented by means of a function G introduced in [8,9]. The 
generalization to a / was achieved in [2]. Using the above notation, we can define a 
generalized G-function and its 'dual' G through the linear integral equations 



(2.22) G = /right + K a -k G , 

(2.23) G = f lc{t + G*K a . 
A formal solution utilizing the resolvent is 

(2-24) G(C, C'K «'; «) = (/right + R * /right) (C, CO , 

(2.25) G(C, C'l«. «'; «) = (/left + /left * R) (C, C) • 



Up to some factor we reproduce the functions G and 

*(C, C\k, a) = (/ left * G) (C, CO = (G* /right) (C, CO • 

from [2] for k' = a + k. In the general case ^ is equal to the terms in brackets on the 
right hand side of (2.8). 

The important 'normalization condition' for the function u, originally suggested 
in [10], includes the asymptotics for large spectral parameter 

(2.26) lim C a Lu((,('\K,K + a;a) = 0. 

In [3] we showed this directly from the definition of the previous section using the 
properties of the linear integral equation (2.22). The key identity in our proof was 

(2.27) lim 2(C/C')- a G(C, CO = Q~ a ~ p(C') ~ I G( V , C')\{vKT a = , 
valid for n' = n + a in the measure. 
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2.4. Functions a and <J>. In order to prove (2.27) we introduced a generalized 
'dressed charge' a and a function $ in [3]. Here we will use slightly modified and 
generalized definitions of these functions adapted to the notation in [5]. Namely, we 
define a function a which satisfies the following linear integral equation 

(2.28) a(C|«, a) = A C C" Q + I dmfa) a( V \K, a)K a (rj/Q . 

The so-called 'dressed function trick' implies that 

/ dm(ry) G(?7, C'|«, oOA^ - " = / dm(ri) (t(t]\k, k'; a)f Tig u(rj, ■ 

Then (2.27) is reduced to an identity for a with k! = k + a which can be proved with 
the help of the properties of the function 

Q(C,k') 



(2.29) *(C|k,k') 



Q(C,«) 



The key observation here is that two TQ-equations (2.18) with twist parameters k 
or k' , respectively, can be combined into the following relation for 

f0 o^ 1 ${v\k,k') *(w|k,k') 



p(t]\k, k')(1 + o(?7,k)) A^$(i?|k,k') p(?7|/«,K')A^$(r7|/«,K') 
Introducing an operator 

(2 - 31) H * = 1 + 4,-) 4 + l + ak«) ^ " ' 

where d^f(rj) = f{q ±1 i]) and a(r),K) = l/a(rj, k), we can rewrite this more compactly 
as 

(2.32) H v $(ti\k,k , ) = 0. 

This is the operator H which we mentioned already in the introduction and which 
will be useful for the characterization of our main function u as well. 

For k! = k + a the solution of the integral equation (2.28) for a can be expressed in 
terms of the function $: 

(2.33) a(() = A c $(C)/<Do 

with some known normalization constant $o [3]- This was shown in [3] for the 
inhomogeneous model with finite N. Here we shall repeat the argument for the 
temperature case. 

For the proof we need to know the location of the singularities of the involved 
functions. If N is finite the function has poles at the Bethe roots corresponding to 
the zeros of the denominator Q(C, «)• It follows from our definition of the Q-operator 
described in Section 3 of [5] that 

(2.34) Q((,k) = C k+s A((,k), 

where A((, k) is a rational function in ( 2 . As explained above we consider the sector 
s = here. Thus, 

(2.35) m^^')/^o = c K ' +K H^^'), 

where A = log£ by definition and where 4>(X\k,k') has constant asymptotics for 
Re A — > ±oo. 
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It is well known that in the temperature case all Bethe roots 'condense' to the point 
C = 1 which becomes an essential singularity of k'). We conclude that (ft(\\n, n') 

has an expansion of the form 

(2.36) 4>(X\k,k') = l + ^ Cj cth^(A), 

i>o 

where cth^(A) := (d/dXy cth(A) and where the dependence on the inverse tempera- 
ture f3 is in the coefficients Cj which also depend on q, k, k'. Note that all partial sums 
in (2.36) are still meromorphic in the original variable £ 2 - The coefficients Cj have the 
high-temperature expansions 

(2.37) Cj = ^ Cfcb . 

k>j 

with respect to j3. The most convenient way to calculate the coefficients c^y is to 
consider the functional equation 

(2>38) $(?CI«,k') _a(C,«') 



$(q-\\K,Ki) o(C,«) ' 

where the high-temperature expansion on the right hand side follows directly form the 
integral equation (2.21). In Appendix C we will explain how to calculate the lowest 
coefficients explicitly. 

For k' = k + a we have &(i]\k, k + a)/<I>o = V~°' ( P( r l\ K i K + Q )> an d the factor rf a 
cancels the corresponding factor under the integral in (2.28) stemming from the kernel 
K a (r)/(). Then, if we substitute (2.30) into the integral and take into account that 
both, the function a(() and the second term in the right hand side of (2.30), do not 
have singularities inside the contour, we conclude that this second term does not 
contribute to the integral. Hence, it can be dropped and the integral can be calculated 
by deforming the contour as it was explained in [3]. 

3. Solving the generalized dressed charge equation by means of the 

high-temperature expansion 

We are now going to generalize the above ideas step by step. We will show that all 
the linear integral equations considered above, like equation (2.28) for the function a, 
equations (2.22), (2.23) for the G-functions or equation (2.15) for the resolvent, can 
be solved in terms of a single function of one or two spectral parameters with essential 
singularities if one of the spectral parameters equals 1 and a certain behavior at and 
oo, as in the above example of the function a(Q at n' = k + a. We shall see in the 
next section that one such function <£>(£, £') with essential singularities at £ = (' = 1 
will be all we need. The full information about the solutions of the above mentioned 
linear integral equations is contained in this single function. 

Before we proceed to the general case let us describe how our method works for 
equation (2.28) with independent parameters k, h' and a. We have seen above that 
a given by (2.33) with 3>(£|k, k') defined by (2.29) solves equation (2.28) only if 
k' = k + a. The problem with the general case is that the behavior at ( = 0, oo 
does not match the corresponding behavior of the kernel K a in the integral of (2.28). 
Hence, we need a generalized function <3?(£|k, k'; a). 

In order for a formula similar to (2.33) to work, we first of all need a relation similar 
to (2.30): 

m ,\ 1 ^( r l\ K i k';o) , . , . 

(3-1) / = ~ a it, i , — s +rvq \k,k ;a) , 

p(ri\ K , K ')(l + a(ri,K)) A„$(»/|/e, a) 
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where the remainder function r is not defined yet. The only property of this function 
is that it does not have singularities inside the contour 7. We will call such functions 
'regular'. We also need a generalization of the asymptotic property and the expansion 
(2.35), (2.36), 

(3.2) $(C|K,K';«) = r a 0(A|K, K';a), A = logC 
with 

(3.3) 0(A|/e, «';<*) = l + J^cth^A), 

where the coefficients c'j now depend on q, k, k', a and on the inverse temperature /3. 
The high-temperature expansion is of the same form as (2.37), 

k>j 

The relations (2.30) and (3.1) can be combined into 

$(i]\k,k') $(r]\K,K';a) , 

with a regular remainder r. Our first claim is that this equation has a unique solution 
in terms of the expansion (3.3). Our second claim is that, due to similar arguments as 
at the end of the previous Section, the solution to the equation (2.28) can be expressed 
as 

(3.6) < T(C|K,K / ;a) = A c *(C|K,/e / ;a). 

We can interpret (3.5) as a cancellation condition for the singular part at A = 
of the functions at the left hand side. Inserting (2.35) and (3.2) we find that the 
following expression is regular (which we denote by '= reg'), 

(3.7) 4>(X\k, n')(q~ a <j)(X + 7rz4|/c, k; a) - q a (f)(X - Trvi\K, k'; a)) 

- 4>{X\k, k; a)(q~ K ' +K (j)(X + -kvi\k, k) - q K '~ K (f>(X - ttvi\k, k')) = reg . 

We substitute the expansions (2.36), (3.3) into here and use the elementary formula 
cth(A) cth(A') = 1 + cth(A) cth(A' - A) - cth(A') cth(A' - A) 

to extract the singular part with respect to A for both terms on the left hand side of 

(3.7) . Comparing the coefficients in front of corresponding terms cth^(A) we obtain 
the following equation for Cj , dj , 

(3.8) (q a - q- a ) Cj - {q K '~ K - q- K ' +K )c'j = 

j 'id /; ^o 

where 

(3.9) 7j (a) := (q a + (-l) j q~ a ) cbh.V\irvi) 

and where the binomial coefficients are defined in such a way that (£) = for a < b. 
If k' = n + a the unique solution is c'j = Cj, since both sides of (3.8) vanish identically. 



11 



Further inserting the high-temperature expansions (2.37), (3.4) we obtain 
(3.10) (q a - q- a )c klj - (q^ - q'^)^ = 

k-l k'-l k-k'-l I" / A / "\ 

E E E c k>\3'Ck-ki\3» {- l ) 3 \\p3'+?'-M)-(- l ) 3 \\p?+r'-M'-K) 

k'=l j '=0 j "=0 I \J / \JJ 

for < j < k - 1. 

If the coefficients Cj are known up to a certain order k of the high-temperature 
expansion (2.37) we can get the c'j from the above relation (3.10) up to the same 
order, because all coefficients on the right hand side are of lower order. Thus, one 
can iteratively solve equation (3.10) with respect to the coefficients c'j, once the Cj 
are known. The solution obtained in this way is obviously unique. In Appendix D we 
show the lowest order coefficients explicitly. 

Remark. Instead of the relation (3.5) we could consider, for example, 

$(r]\K,K';a) _ k'; -a) = 

A v ^(ri\K,K';a) A^$(r?|K, k'; -a) 

Then, repeating the above arguments, we come to the following compatibility condition 

(3.12) (g«_g-«)( C /.+^.) = - £ (-l) J 'X4-''((- 1 ) J t{V(0) 7 ^'''- J ' (O) ^ 

j i j /'>o ^ \ j / \ J / / 

where the coefficients c'j appear in 

(3.13) <t>(\\K,K';-a) = 1 + J^gJ- cth^(A) 

and have the high-temperature expansion 

(3-14) ^ = E^ b - 

k>j 

4. Resolvent and master function 

We shall show in this section that the arguments applied above in order to obtain a 
solution a of the dressed charge equation (2.28) in terms of a function which has a 
simple high-temperature expansion, can be generalized to the linear integral equation 
(2.15) for the resolvent R. To this end we introduce another function <&(£, C'I K > K 'j a ) 
now depending on two spectral parameters (, As we shall see, this function contains 
the whole information which is necessary for the description of all functions defined 
by means of linear integral equations that appeared before. For this reason it may be 
called a 'master function'. 

Generalizing the ideas of the previous section we define it as follows: 

(4.1) $(c, Ok, «) = Cl«, «) + a^V(C/C, a) , 

(4.2) $ , (C,Cl«,«';a) = i(C/C)XA,A , |«,«';a), A = log (C) , A' = log(C), 
where 

(4.3) ^(A,A'|K,*e';a) = c^v cth^(A) cth^(A') 
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and where A^ip^, a) was defined in (2.13). We complete the definition by demanding 
that <]?(£, C'I K > K '; a ) satisfies an equation similar to (3.1), 

u 1 _ &(v,('\k,k';<x) ( r i\ i \ 

l4 - 4j p(r,\ K , k')(1 + afa, «)) " A^fa, C'|«, tf; a) + rW ' " j 

where the remainders r, r' are regular functions of 77. 

We shall see below that (4.2) and (4.4), (4.5) determine the coefficients Cjj> in 
the high-temperature expansion (4.3). The relation of our master function with the 
resolvent is explained in the following 

Lemma 4.1. The resolvent R defined by the linear integral equation (2.15) can be 
represented as 

(4.6) R((, C'\k, a) = - J-,A C A C <D(C, CV, «'; <*) • 
A proof of this Lemma is provided in Appendix A. 

The function <3? of one spectral parameter introduced in the previous Section is 
recovered in the limits of large (, 

(4.7) lim C a R((, C'\k, «'; a) = --^(Cl*, «'; a) = _-La ( ,*(C'|k, a) , 

C^oo 47T7 47T7 

lim C' a R(C, C'\k, a) = --^ct(C|k, k'; -a) = --^A c $(C|k, «'; -a) , 



or 



$(C'|«> a ) 



(4.8) lim C" a *(C,C'|K.«';«) 



lim ("*$((, C'\k, k: a) 

Let us now come back to our claim that the high-temperature expansion (4.3) is 
determined by (4.4), (4.5) and by the asymptotic condition (4.2). Combining (4.4), 
(4.5) with (3.1) we obtain 

(A q\ $>(ri\K,K';-a) $'(77, ('\k, k'; a) _ , , , 

1 ' A v $(ri\K,K>;-a) A„$(r/, ('\k, k>; a) r W^I K '"'^' 

, A -,M $>(r)\K,K';a) $'(£, i]\k, k'; a) ~ , 

(4-10) - — , ; r — ; : ; r = T 1 (C , 77 K, K J a) , 

A„$(»/|K,/e';a) A„$(C,»7|K,K';a) 

where the remainders r, f' are also regular with respect to the variable n. 

First of all we shall see that (4.9) and (4.10) are equivalent to each other. We 
substitute (4.1) into (4.10) and obtain the regularity condition 

(4.11) 

4>(h\k, k'; a) ^q~ a (f>(\, fi + TTvi\n, k'; a) — q a (f>(X, \i — irvi\K, k'; a) — cth(A — 

— <p{\, fi\K, k; a) ^q~ a 4>([i + TTvi\n, k; a) — q a (f)([i — nui\K, k; a)^j = reg , 

where /j, = log (77) and where the variable £ corresponding to A = log (() must be 
outside the integration contour 7. 

If we repeat the procedure described in the previous Section, namely, if we substitute 
the expansions (3.3), (4.3) into (4.11) and into the equation obtained from (4.11) by 
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replacing a — > — a and then take the singular part, we obtain two equations for the 
expansion coefficients Cjj/ 

(4.12) c ^'Uf,l = 
j'>o 

(4.13) UjA-iy'^fj = -(-1)' 
j'>0 

where 

(4-14) U j , l = (q a -q- a )5 j>l + ^ > 

j'>0 

(4-15) U jJt = (q a - q~ a )kl ~ E (-^"'^'b 

J'>0 

and 

(4-16) 7^11 = ((" 1 ) J (f) " (-l/Q) 

with 7j defined in (3.9). 

The compatibility condition for the equations (4.12), (4.13) is 

i'>o \ 3 / r > V ( / 

Substituting (4.14) and (4.15) here, using the identity 
(4-17) 

where 

(4-18) 7i'j»|i = (({) " (-l)'(f)) 

and doing some algebra, we come back to our previous compatibility condition (3.12). 

Hence, either of the two equations (4.12) or (4.13) can be used to determine the 
coefficients Cjj>. Their high-temperature expansion is of the form 

(4.19) c jd ,= J2 P kc k\j,j>- 

k>j+j' 

Substituting this together with (3.4), (4.14) into equation (4.12) we obtain 

/ ■ _|A k-i k'-j-ik-k'-i 

(4.20) (q a -q~ a )c m = ( J ^ £ c k'\j,j' c k-k'\j"^',j"\l ■ 

V J k'=j+l j'=0 j"=0 

Again, like in equation (3.10), the coefficients c^uy are completely determined by the 
previous coefficients c k i^/ y with k! < k. And again we obtain an iterative calculation 
scheme which is very efficient and allows us to quickly calculate high-order terms in 
the high-temperature expansion (4.19). 
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5. Representation of G, G and uj in terms of the master function 

Since we have the formulae (2.24), (2.25) and (2.8) which relate the functions G, G 
and uj, respectively, to the resolvent R given by (4.6), it is merely a technical problem 
to express these three functions in terms of <!>(£, £'| k j n'; a). An important point here 
is that it turns out to be possible to get rid of all integrations. 

Lemma 5.1. The functions G, G and uj are expressed in terms of the function <3? as 
follows 

(5.1) G((, C'k «'; a) = A c fl- c ^(C, C'k «'; a) , 

(5.2) G(C, C'\k, a) = - 2^fT c A c ^(C, C'k «'; «) , 

(5.3) ^(C, C'K «'; a) = H C H C $((, C\k, k'; a) , 

where the operators A^ and are defined in (2.10) and (2.31), respectively. 

The proof of this Lemma is deferred to Appendix B. 
The master functions $ has the symmetry 

(5.4) <!>((,('\ K ,K',a) = <!>((',(\K, K ',-a) 

which follows from (4.1), (4.9) and (4.10). Due to our Lemmas 4.1, 5.1 this symmetry 
carries over to R, uj, G and G, 

(5.5) R((, C'\ K i K> ') a ) = -R(C') CI K > K '> ~ a ) > 

uj{(, C\k, k', a) = uj{C, (\k, k', -a) , 

G(C, C'|«, a) = " 2^G(C', Cl«, "«) • 

A few comments are in order here. First, the expression (5.3) for our main function 
uj is rather appealing as compared to the original expression (2.8). We could get rid 
of all unpleasant integrations that make the analysis hard. Second, we derived our 
formulae in the temperature case. Yet, (4.6), (5.1), (5.2) and (5.3) are valid in the 
general inhomogeneous case for finite N. This means that we can use our formulae 
for the future analysis of the CFT scaling limit. 

As we pointed out above our master function $ of two spectral parameters encodes 
the complete information about various functions appearing in the description of the 
physical part of the correlation functions of the XXZ chain. This function $ can be 
characterized by either of the relations (4.4), (4.5), (4.9) or (4.10) which ensure the 
cancellation of singularities inside the integration contour 7. These relations make 
sense in the general case as well, when the lattice is finite in Matsubara direction. In 
that case they mean that the residues at all poles corresponding to the Bethe roots 
must vanish. In the temperature case we could use the high-temperature expansion 
(4.3) and determine the coefficients Cjji. This procedure can be efficiently implemented 
on a computer. We hope that it will be possible in the future to find a generalization 
of this procedure to the CFT scaling limit, where relations like (4.4), (4.5) or (4.9), 
(4.10) may turn into a Riemann-Hilbert problem. 

We also believe that it would be interesting to understand the deeper meaning of 
the operator H defined in (2.31). 

Acknowledgments. The authors are indebted to M. Jimbo, Kh. Nirov and F. Smirnov 
for helpful discussions. HB is grateful to the Volkswagen Foundation for financial 
support. 
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Appendix A. Proof of Lemma 4.1 

The statement of the Lemma can be verified rather directly. First we write (2.15) 
in the following way, 

(A.l) R(C, = K a ((/C) + (K a * R)(C, CO , 

where we omit the other arguments of the resolvent for simplicity. We will also omit 
such arguments in all other functions during this proof. 

Let us substitute the right hand side of (4.6) into the right hand side of (A.l) and 
use (4.4). Then 

J 1 T p{v)(l + a(r])) 

= KM/C) + jf ^«.»>(^7) - >-fer))^A,(#(„, 9 c') - 

where we used the arbitrariness of the parameter r and choose either r = q(' or 
r = Q^C- We further rewrite the last line as 

K a {C,K') + J ^r*«(C, V)^-^ (*'(r,, C) - r(r,, C)A„<l>fa, Co) 
= K a (C/C') + ^ jf ^«(C, ^A^fa, CO 

= *«(C/C) - ^A c A c d>'(C,C0 = ^A c v(C/C',«) - ^A c A c ,d>'(C,C0 
= -^(^(CO + V(C/C») = -^A f A^ (C ,CO , 

where we deformed the contour in such a way that it does not contain the essential 
singularity rj = 1, but the two poles coming from the kernel K a (^,r]) instead. One 
also has to verify that the contribution of the boundary terms is equal to 0. This can 
be done in a similar way as it was explained in [3]. 

Appendix B. Proof of Lemma 5.1 

(i) Let us start with (5.1). In accordance with the formula (2.24) we have (again 
we keep only spectral parameters in arguments and omit all other arguments) 

G((, CO = /right(C, CO + (R * /rightXC, CO , 

/right (C, CO = v>(rt/cO - p(('M(/0 ■ 

Hence, the integral on the right hand side is 
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where both spectral parameters ( and £' are inside the integration contour 7 and 
[Arf RK,V) ,, -1 . gt^j) ,, 

h= j 1 ^7m^m* lq " /c) ' I ' = J^7m^m* l ' ,K) - 

We treat the integral ii in a similar way as in the proof of Lemma 4.1, 

- A,$h-'C,»)(- + ' , [f'CD))i'fa-VO ■ 

The integrals containing the remainder function r' do not contribute, because all 
singularities are outside the integration contour. We conclude that 

drj 2 1 



rj 2 2ni 



&{qC,v) - ^V'C, VO = A c d>'(C,<zC0 



In the second integral I2 we have to take into account the residue originating from 

/(')■■ 

h = -%o(i + «(o) + 1 iFim+m*^ ] ' 

where the contour 7' does not contain the point £'• The corresponding integral can be 
calculated in the same manner as I\. The result is 

Altogether 

g(cc') = y^c/co - Picm/o +h- P (Ci)h 

= iKrt/O - piCm/C) + a c $'(C, qC) + 2m Y^r) - ^')A c <f'(C, CO 
= a c $(c, gO - p(OAc*(C, CO - T - 7 -^yA c A c $(C, CO 

(ii) The proof of the formula (5.2) proceeds similarly to the above proof of (5.1). 
We end up with 

B(CC') = - ^a^a^k.O = -±. Hi A ( , m o. 

(iii) Finally we prove (5.3) which is the main result of this paper. The function cu is 
given by (2.8) 

(B.l) ^(CCO = (G* /right)(C,CO + 5 c -<5 f -A c -V(C/0 • 
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After inserting (5.2) the first term on the right hand takes the form 

f An 2 1 1 
(R2) " / ^-^ p(,)(l + a(,)) W(C,^(i/, O = 

1 /(+) 1 /(-) | o(C)I iO) | g(Q , P(C0 7 (-) m ^UO) 

l + a(C) 71 l + a(C) 71 +PlC)/l + l + a(c/ 2 + l + a(C) 2 P(CMC } ' 2 ' 



where 



In order to calculate all these integrals we proceed in a similar way as before 



Since £' is inside the contour one can see that the first term inside the bracket does 
not contribute, and we have 

(R4> = hi / ^k^'c-jja^c^'cjvwo 

= - r( ,«c, owe. o = -*V'c, o - ^Cct ' 

Furthermore, 

= -^ , (C,9C)-KC,C)V'(?" 1 C/0 



18 HERMANN BOOS 1 "* AND FRANK GOHMANNt 

since the term with lim£»_^ vanishes. Finally 

= -r(c,C / )A c ^(C,O-KC,0V'(C/C / ) 

= -$'( C co - Ac,(I)(c ' c/) 1 <KC/C) 

^ j p(C')(l + a(C0) p(C)(l + a(C))^^ /U ' 

If we now substitute (B.3), (B.4), (B.5), (B.6) into the expression (B.2), combine 
with the last term at the right hand side of (B.l) and do some algebra, we come to 
the final result (5.3). 



Appendix C. Coefficients in the high-temperature expansion of the 

FUNCTION 4>(X\k, k') 

As we mentioned in Section 2.4, we can compute the coefficients Cj of the expansion 
(2.36) by solving the functional equation (2.38) and by using the corresponding high- 
temperature expansion (2.37). Unfortunately, the resulting formulae for the coefficients 
Cfc|j are not simple, not even at lower orders. There is a better object though, for 
which the coefficients of the high-temperature expansion are much simpler and from 
which the coefficients cuj can be obtained. 

It might seem that the function A considered in Section 2.4 is a candidate to be 
such simple function. Unfortunately, A is not well-defined in the temperature case 
when the limit N — > oo is performed. Therefore we need another function A which 
exists in this limit. In order to define it we first introduce a function (f) by changing 
the normalization of (j), 

(C.l) $(C|K, K / ) = r K ' + >(A| K ,K'), A = logC. 

Due to (2.36) the function <f> has an expansion of the form 

(C.2) 0(A|«,«') = l + X)c j cthO")(A), c j = Y i Pc k \ j . 

j>0 k>j 

It can be obtained from the expansion (2.36) by taking into account that $o = 
(4>(cc\k, k') +</>(— oo\k, k'))/2. Then there is a unique solution A(X, k) of the functional 
equation 

(C.3) 0(A|k,«') - A(A ' K,) 



A(\,k) 

in form of a high-temperature series 

(C.4) 2(A, k) = expi ^ Oj coth( J ')(A) I , aj = ^ P k a k]j , 

i>o ' k>j 



with coefficients a k y which depend only on n. These coefficients are comparatively 
simple. 



Here are the first few of them, 



ft T 1 /"I 

u l|0 


-fa - 


9 7/ ^5 


?; — a 2n + 1 




a 2|0 = 


-(a 2 - 
19 


9 )\ u 






0-211 


-(a - 


9 7 V 17 






a 3|0 = 




9 A w 


l)(t;-2)((g 2 + g- 2 )( W 2 - 


- 6u + 6) + Av 2 - 12v + 12) /(3t> 5 ) 


"311 


-(a - 


9 7l9 


g-2 )(w _ 1)(t; _ 2) 2 /(27; 4 : 


), 


a 3|2 = 


-in - 


9 J 


l)(t;-2)/(12^), 




«4|0 = 




-9' 2 )(^- 


l)(t;-2)((g 2 + g- 2 )( W 4 - 


- I8v 3 + 78w 2 - 120w + 60) 








+ 8(^ 4 


- 9v 3 + 24t> 2 - 30v + 15)) /(3v 7 ) , 


a 4 |! = 


-(9- 




1)((9 2 + 9 _2 )(11^ 4 -138 


v 3 + 498v 2 - 720v + 360) 








+ 2(l9v 4 - I86v 3 + 546t> 2 - 720v + 360)) /(36w 6 ) , 


«4|2 = 


-(9- 


9- 1 )V " 


- q~ 2 )(v - l)(v - 2){v 2 - 


6w + 6)/(12t; 5 ), 


«4|3 = 


-(9- 


9- 1 ) 4 (^- 


l){v 2 -6w + 6)/(144v 4 ). 





The explicit coefficients c^y can now be obtained following the above prescription 
in opposite direction. 



Appendix D. High-temperature expansion of the function $(C|«, n';a) 

Following the above scheme, one can obtain the coefficients Cj up to any order in (3. 
In this sense we shall imply now that all these coefficients are already known. The 
next step is to calculate the coefficients c'- from (3.3) which determine the function 
3>(£|re, «/; a) given by (3.2). To this end we can simply solve equation (3.8) with 
respect to the c' k ^. Here are the first coefficients obtained this way 

_ _ q a - g- a 

C 1|0/A - Cl| , X ■- ^ K ,_ K _ ^_ K ' +K i 
C 2|o/X = c 2|0 

( ? + ( g («'-«+")/2 _ g -(«/-K+a)/2) ^(k'-k-«)/2 _ ? -( K '- K -a)/2) ^ 

+ ( - g -i)( «'-«- g -«'+«) Cl| ° ' 

4|l/X = C 2|l > 
c 3,o/X = C 3 |0 

(q + g -l)(g(«'-«+")/2 _ q -(it'-it+a)/2j ^(«'-«-a)/2 _ ? -( K '-«- a )/2^ 

+ 2 (9-9- 1 )(9 K '- K -9- K ' +K ) Cl| ° C2| °' 

( g + g-l)2(q(K'-K+a)/2 _ g -(«'-K+a)/2^2^( K '_«_ Q )/2 _ ^-^'-^-a)^ ^ 

+ (9-9- 1 ) 2 (9 K '- K -9~ K ' +K ) 2 Cl| °' 

C 3|l/X = C 3 |l 



( g + g -l)( 9 (K'-K+a)/2 _ q -' K '- K+a )/2j ^ q ( K '- K - a )/2 _ q -( K > - K - a )/2j 
C 3\2/X = C 3 |2- 
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Appendix E. High-temperature coefficients for ('\k, k'; a) 

Here we show some coefficients Cjj>,c k \j ji that determine the function ('\k, k'; a) 
via the expansions (4.3) and (4.19). Let us start with two useful observations: 

Remark E.l. Analyzing the equations (4.12), (4.13), we observe the following 

c' d 
(E.l) Cj o = , c i = . 

v ' J' u qa g— a ' v " q— a 

This can also be deduced from the limiting relations (4.8). As in the Remark above, 
these two relations hold to all orders in the high-temperature expansions. 

Remark E.2. The equations (4.12), (4.13) have the symmetry 

(E.2) c ;. _ ( _i)^, 

Cj j , (_iy+3' Cjl j . 

It is easy to check that the equation (3.12) is also explicitly symmetric under the 
transformation (E.2). 

Due to Remark E.l it is easy to obtain the high-temperature expansions for Cjo 
and co ,j. We will not show the corresponding coefficients c^q and c k \o,j here. The 
other coefficients can be obtained solving equation (4.20) order by order. Here are 
several examples: 

2 , 

^l 1 . 1 _ q a _ q-a C 3|2> 

_2 , 2{q + q' 1 ){q a + q' a ) , , (q + q" 1 )^ + q~ a ) , 2 

C 4|l-1 " q a _ q-a C 4|2 + (? _ q - 1){g a _ q - a)2 C 1|0 C 3|2 {g _ _ Q - a)2 ^ C 2\l) , 

3 , 

C4|l,2 - C4|2,l - ga _ c 4|3- 

All the coefficients c^,, . can be obtained following the previous Appendix. With the 
help of a computer one can efficiently obtain higher coefficients c' k ^ and c k yj/ as well 
by means of the above scheme. 
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